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(4,4')-Bipyridine in vacuo and in solvents: A quantum chemical study 
of a prototypical floppy molecule from a molecular transport perspec¬ 
tive 
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We report results of quantum chemical calculations for the neutral and anionic species of (4,4')-bipyridine (44BPY), a proto¬ 
typical molecule with a floppy degree of freedom, placed in vacuo and in solvents. In addition to equilibrium geometries and 
vibrational frequencies and spectra, we present adiabatic energy curves for the vibrational modes with significant intramolecular 
reorganization upon charge transfer. Special attention is paid to the floppy strongly anharmonic degree of freedom of 44BPY, 
which is related to the most salient structural feature, namely the twist angle 9 between the two pyridine rings. The relevance 
of the present results for molecular transport will be emphasized. We show that the solvent acts as a selective gate electrode and 
propose a scissor operator to account for solvent effects on molecular transport. Our result on the conductance G vs. cos^ 9 is 
consistent with a significant transmission in perpendicular conformation indicated by previous microscopic analysis. 


1 Introduction 

The present work represents a quantum chemical study in¬ 
tended to pave the way toward microscopic nanotransport 
studies in single-molecule devices in vacuo and in solvents 
based on (4, 4')-bipyridine (44BPY), a prototypical floppy 
molecule, wherein a floppy vibrational mode turns out to yield 
an important intramolecular reorganization upon charge trans¬ 
fer. This is an effect, which existing transport studies did not 
consider. 

In view of its special structure, 44BPY is a representa¬ 
tive model molecule that attracted considerable experimental 
and theoretical interest. Compounds based on 44BPY, com¬ 
monly known as viologens, captured chemists’ attention for 
many decades, especially due to their electrochemically re¬ 
versible behavior.^ They have been investigated from various 
perspectives, like herbicidal activity,— electrochemical display 
devices, storage energy, catalytic oxidation,—"— photoactive 
molecular, and supramolecular machines.— The preparation of 
complexes embedding neutral, radical anionic and dianionic 
forms of (4, 4')-bipyridyl ligands is of interest for design¬ 
ing switchable materials wherein the redox chemistry can be 
finely matched to facilitate electron transfer. 

More recently, the 44BPY molecule has been used in 
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molecular electronics.^— Due to its special structure, with 
two active nitrogen atoms in para positions (cf. Fig.[T]), 44BPY 
is particularly suitable for simultaneous binding to metallic 
electrodes. This fact has been exploited in the first work that 
succeeded in the repeated formation of numerous molecular 
junctions.— 44BPY molecules have been incorporated in re¬ 
dox active tunneling junctions to demonstrate the electrolyte 
gating.— 

It is from the perspective of molecular electron transport 
that the present quantum chemical study has been conducted. 
Unlike in most molecular electronic devices fabricated so far, 
in the single-molecule junctions based on 44BPY there exists 
clear experimental evidence that the electric current is due to 
negatively charged carriers (n-type conduction); the measured 
Seebeck coefficient is negative.— This indicates that the elec¬ 
tron tunneling through molecular junctions wherein a 44BPY 
molecule is contacted to gold electrodes is predominantly due 
to the lowest unoccupied molecular orbital (LUMO). A posi¬ 
tive Seebeck coefficient, with the associated p-type conduc¬ 
tion, would have revealed the predominance of the highest 
occupied molecular orbital (HOMO). As a further support, 
one can invoke the results on electrochemical gate-controlled 
electron transport through redox-active viologen molecular 
junctions,— which also indicate a LUMO-mediated conduc¬ 
tion. This kind of experiments also enables one to distinguish 
between electron (LUMO-mediated) n-type conduction and 
hole (HOMO-mediated) p-type conduction: the gate poten¬ 
tial (overpotential in electrochemists’ nomenclature) shifts the 
energies of the occupied and unoccupied molecular orbital in 
opposite directions relative to the electrodes’ Fermi level, en¬ 
hancing or diminishing the electric current. 





By implication, this means that, when an electron subjected 
to a source-drain voltage travels across a 44BPY-based molec¬ 
ular junction, a transient radical anion (44BPY*^) is formed. 
Therefore, the determination of the properties of the neutral 
44BPY‘’ molecule as well as of the bipyridyl radical 44BPY*^ 
is a prerequisite for understanding the electron transport in 
these molecular junctions. This is the aim of the present paper. 


The quantum chemical calculations reported below have 
been performed by using the Gaussian 09 suite of programsi^ 
at DFT/B3LYP level. Our study substantially goes beyond 
earlier works wherein both the parent molecule TTBPY*' and 
the radical anion 44BPY*^ have been investigated theoreti¬ 
cally, the latter especially in relation with photoreduction.i^— 
The aforementioned theoretical works used smaller basis sets 
and only considered a few ground state properties of 44BPY*^ 
and 44BPY*^ in vacuo. In addition to the case of vacuum, 
we also present results for 44BPY in solvents of experimental 
interest. The solvent has been described within the polarized 
continuum model using the integral equation formalism (key¬ 
word SCRF=IEFPCM in Gaussian 09). For the anionic (and 
cationic, cf. Sec. |5]l species, we carried out unrestricted cal¬ 
culations (Gaussian keyword UB3FYP). Because calculations 
for open shell systems can sometimes be problematic (wrong 
geometries, spurious symmetry breaking, or spin contamina¬ 
tion, see, e. g.. Ref. l2ll) . to avoid artefacts we always checked 
our results carefully, and convinced ourselves that, for in¬ 
stance, optimizations properly converged (all frequencies are 
real), symmetries (point group D 2 for44BPY’’ and 44BPY*+, 
and D 2 /J for 44BPY*^) and spin are correct. For concreteness, 
we mention that, with the default settings of Gaussian 09, typ¬ 
ical spin values found are {S^) = 0.7592 before annihilation 
of the first contaminant and (5^) = 0.7501 after annihilation. 
To be complete, let us still note that no problem arose even in 
the “small” calculation for the anion done at the very modest 
ROHF/6-31g(H-*) level of theory (cf. end of Sec. [2i; that cal¬ 
culation has been solely done to allow comparison with and 
check results of earlier work. 


The remaining part of this paper is organized as follows. In 
SectionlH results for the equilibrium geometries of the neutral 
and anionic species obtained by using basis sets larger than 
those employed in previous studies will be reported. Further, 
vibrational effects in the radical anion placed in vacuo and in 
experimental conditions (acetonitrile as solvent—) will be dis¬ 
cussed in Section|3] Section|4]is devoted to an analysis of the 
relevant adiabatic energy curves for 44BPY*^ and 44BPY*^; 
this analysis is augmented with the discussion in the Elec¬ 
tronic Supplementary Information. Results for the 44BPY 
molecule immersed in solvents will be presented in Section 
|5] Because the present quantum chemical investigation is in¬ 
tended to pave the way toward transport studies in 44BPY- 
based molecular devices, we present in Section|6]preliminary 
results revealing the impact of the molecular conformation on 


the electrical conductance. Discussions and conclusions will 
be presented in Section|2l 

2 Equilibrium geomerties of 44BPY® and 
44BPY* 

In this study, quantum chemical calculations for geometry op¬ 
timizations without imposing any constraints have been per¬ 
formed both for the neutral and for the anionic species. In such 
calculations, an important issue is to employ a sufficiently 
large basis set. 

We have extensively investigated the impact of the basis set 
size. Table SI (tables with label S are from the Electronic 
Supplementary Information) collects our results on the equi¬ 
librium geometry of the neutral 44BPY molecule obtained 
via DFT/B3FYP and various basis sets. The atom number¬ 
ing is given in Fig. [T] The optimum geometry of the neutral 



Fig. 1 The (4,4') molecule (44BPY) and the atom numbering used 
in the present work. 

species has been previously investigated in a series of stud- 
ies. 13,14,17,18,20,22 xjigse studies used lower theoretical lev¬ 
els (semi-empirical, HF) and/or smaller basis sets (sto-3g, 3- 
21g(H-*), 6-31g(H-*)^). As compared to those results, our re¬ 
sults represent an improved agreement with experiment .212^ 
For the cc-pVTZ basis set, the largest set used earlier for 
44BPY‘',— our study basically confirms the results reported in 
Ref. [ 1 ^ We note that there exist small differences between the 
present geometrical parameter values computed with Gaus¬ 
sian 09 and those (also included in Table SI) computed with 
GAMESS in Ref. Hat the same DFT/B3FYP/cc-pVTZ level 











of theory. However, they are small; bond lengths differ by at 
most 0.001 A. By and large, Table SI confirms the previous 
claim— that the geometrical parameters of 44BPY‘’ are little 
sensitive to the basis set size. An important fact to be em¬ 
phasized in connection to Table S1 is that, while substantially 
increasing the computational cost, the usage of larger basis 
sets does not necessarily improve the agreement between the¬ 
oretical and experimental bond metric data. For instance, al¬ 
though the inter-ring C-C bond length is somewhat reduced 
by increasing the basis set, it remains slightly larger than the 
experimental one.— The theoretical value becomes closer to a 
normal C(sp^)-C(sp^) single bond (1.479 A).— 

Most sensitive to the basis set size is the twisting angle 
9 = C 3 C 1 C 2 C 4 between the two pyridyl rings, which has been 
long recognized to be the salient structural feature of this 
molecule, 13,14,17,18 ^vjjose variations are related to very small 
energy variations. 

Our bond metric data for the radical anion 44BPY*^ are 
shown in Table S2. One should note that the previous the¬ 
oretical studies on the anion have been done either semi- 
empirically— or by using much smaller basis sets (HF/3- 
21g(H-*)iI and DFT/6-31g(H-*).i^ In the absence of direct 
experimental results for the anion structure, we can only 
present the comparison with existing theoretical results ob¬ 
tained for the largest basis, namely 6-31 g(H-*) of Ref. [18. 
Somewhat surprisingly, the geometrical parameters are little 
sensitive to the basis set size even in the case of the radical 
anion 44BPY*^. Noteworthy, the computed geometrical pa¬ 
rameters are rather insensitive to the presence of the diffuse 
basis functions. This behavior contrasts with that exhibited by 
other properties of the anion, which will be examined later in 
Sections II] and| 6 ] 

Because the structural differences between the neutral 
molecule 44BPY‘' and its reduced form 44BPY*^ have been 
amply discussed in the literature (see, e. g.. Refs. FTlflS and 
m, we only briefly mention the main aspects. The twisting 
angle 0 between the two pyridyl rings results from the inter¬ 
play between the TT-electronic interaction of the pyridyl frag¬ 
ments and the steric repulsion of the ortho hydrogen atoms 
(Hi 3 and H 14 , and H 19 and H 20 in Fig.[T]). The former favors 
a planar conformation ( 0 = 0 ) wherein the ;r-electrons are de¬ 
localized between the pyridyl rings, while the latter prefers a 
nonplanar conformation ( 0 ^ 0 ) wherein the steric hindrance 
is reduced. The neutral molecule can be characterized as con¬ 
sisting of two aromatic pyridyl rings linked by a (nearly) sin¬ 
gle C-C bond in a twisted conformation (point group D 2 ). The 
balance between the two aforementioned effects in the neutral 
molecule yields an equilibrium value Oeq = 37.2°, as measured 
in electron diffraction experiments,— which can be well repro¬ 
duced by quantum chemical calculations (cf. Table SI). 

As suggested for the first time on the basis of a semi- 
empirical MO approach,— and later within HF/3-21g(H-*)i^ 


and DFT/6-31g(H-*)ii^ calculations, two important structural 
changes occur on going from the parent molecule TTBPY** 
to the radical anion 44BPY*^. First, the TT-electronic in¬ 
teraction is enhanced by the extra electron delocalized over 
the reduced molecule, which becomes planar {Oeq = 0 , point 
group D 2 /i). Second, the anion is characterized by a (planar) 
quinoidal structure. Basically, the quinoidal distortion is lo¬ 
calized in the inter-ring region. On going from the neutral to 
the reduced species, mostly affected are the inter-ring C 1 -C 2 
bond length, which shortens by ^ 0.052 A, becoming thereby 
substantially stronger and acquiring partial double bond char¬ 
acter (bond order 1.247, see below), and the neighboring Ci- 
C 3 bond length, which increases by ^ 0.034 A. Less affected 
are the C 5 -N 7 bond, which lengthens by ~ 0.018 A, and the 
C 3 -C 5 ring bond, which shortens by ^ 0.012 A. The calcu¬ 
lated bond orders are only slightly dependent on the method. 
With ROHF/6-31g(H-*) we found 1.288 and 1.532 for C 1 -C 2 - 
and C 3 -C 5 -(Wiberg) bond order indices, confirming thereby 
the previously reported values 1.29 and 1.53, respectively;— 
B3LYP/aug-cc-pVTZ yields the values 1.247 and 1.560, re¬ 
spectively.— 

3 Vibrational modes of the radical anion 
44BPY* 

For the neutral 44BPY molecule, theoretical results obtained 
by using DFT/B3LYP and good basis sets (e. g., cc-pVTZi^) 
have already been compared with the experimental data avail¬ 
able for the Raman— and infra-red (IR) activ e vibrational 
modes. By using basis sets larger (cf. Table SI) than that em¬ 
ployed in Ref. [ijl we did not And a significant improvement 
against experiments; the empirical factors, which are com¬ 
monly introducedto scale the calculated vibrational fre¬ 
quencies, remain close to the value ^ 0.98 found in Ref. [l^ 
Therefore, we confine ourselves to discuss vibrational ef¬ 
fects for the radical anion, wherein earlier resultsiii^ have 
been obtained at modest theoretical levels (HF/3-21g(H-*)— 
and DFT/B3LYP/6-31G(H-*).i^) The results for the frequen¬ 
cies of the Raman active modes, which can be compared with 
available experimental data,— are collected in Table S3. Al¬ 
though Raman measurements have also been carried out for 
the partially or totally deuterated radical anion—), we confine 
ourselves to the undeuterated case, of more direct relevance 
for electron transport. The Raman measurements^ comprise 
the spectral range 300 cm^^ < (0< 2000 cm^*. It does not in¬ 
clude the spectral range (O > 3000 cm^' wherein C-H stretch¬ 
ing modes are active (see Figs. |2] and |3ll. Out of the eight 
Raman active modes (Ag symmetry) shown in Table S3, only 
seven have been detected experimentally.— The lowest ring 
deformation mode ( 6 a, cf. Table S3) has not been observed. 
This can be understood: the spectra calculated by us (Fig. Ell 







show indeed a very weak spectral intensity for this mode. 

As compared with the approach HF/3-21g(+*) of Ref. [B 
where the empirical scaling factor amounts 0.892, the ap¬ 
proach of Ref. [is] already represented a substantial improve¬ 
ment: the empirical scaling factor needed there was substan¬ 
tially larger {fsc = 0.968). The agreement is further improved 
by using larger basis sets, as expressed by the values fsc —0.98 
of Table S3. So, the agreement is basically as good as that for 
the neutral molecule. To conclude, the DFT approach is able 
to accurately reproduce experimental vibrational frequencies: 
the deviations amount to ^ 2%; for the modes shown in Table 
S3 this represents an absolute deviation of ^ 10 —30cm^'. 
One should still note that, similar to the case of the geo¬ 
metrical parameters, further improvement cannot be achieved 
by employing larger basis sets (for aug-cc-pVQZ, we found 
fsc = 0.980), a fact which reflects the intrinsic limitation of 
the current theory. 

In contrast to the fact, already noted, of relatively insensi¬ 
tive vibrational frequencies, the Raman intensities are much 
more sensitive to the basis set. This is exemplified by the Ra¬ 
man spectrum presented in Fig.|2]— 

Because the Raman measurements refer to the anion in sol¬ 
vent (acetonitrile),— we have also performed calculations for 
this situation. As visible in Table S3 and Fig.|2] while weakly 
affecting the vibrational frequencies, the solvent has a signifi¬ 
cant impact on the spectral Raman intensities. 

Although we are not aware of experimental studies on vi¬ 
brational spectra of the infra-red (IR) active modes in the radi¬ 
cal anion 44BPY*^, for completeness we have also calculated 
the IR spectrum. These results are presented in Fig. |3] As 
compared to the Raman spectra, the IR spectra depicted in 
Fig. [3] turn out to be much less affected both by the basis set 
size and by the presence of the solvent. 
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4 Adiabatic ground state energy curves of 
44BPY0 and 44BPY* 

The coherent electron transport through single-molecule junc¬ 
tions based on 44BPY was investigated theoretically in several 
works. In spite of significant differences between them, 

there is a common finding in those studies, namely that the co¬ 
herent electron transport is determined by the LUMO: close to 
the metallic Fermi energy, the calculated transmission exhibits 
a single pronounced Lorentzian peak and the Seebeck coeffi¬ 
cient is negative. The LUMO energy can be obtained from 
electronic structure A-DFT calculations: it represents the dif¬ 
ference between the anion and neutral ground state energies 
Elumo = Ea— Eai at the neutral equilibrium geometry Qn-— 
Refs. FnllB JB assumed a molecular geometry frozen at 
Q = Qn{= 0)- Our recent study— drew attention to the fact 
that intramolecular reorganization is important and should 


Fig. 2 Theoretical Raman spectra calculated with the basis sets 
6-31g(-l-*)^ and aug-cc-pVTZ for the radical anion 44BPY^* 
placed in vacuo and acetonitrile. They illustrate that both the basis 
set size and the solvent significantly affect the Raman spectral 
intensities. The computed spectral lines have been convoluted with 
Lorentzian functions of half-width 20 cm^ *. See the main text for 
details. 


be considered, especially because the 44BPY molecule pos¬ 
sesses a floppy degree of freedom. In particular, this implies 
that the LUMO energy is not fixed at the value Ehjmo = 
Ea{Qn) — E[^{Qm)- Instead, one should consider a distri¬ 
bution of values Elijmo{Q) = £^a(Q) —and perform 
a Q-ensemble averaging of the current computed at various 
molecular geometries Q. The weight function needed for this 
ensemble averaging requires the knowledge of the adiabatic 
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Fig. 3 Theoretical IR spectra calculated with the basis sets 
6-31g(+*)i^ and aug-cc-pVTZ for the radical anion 44BPY^* 
placed in vacuo and acetonitrile. They illustrate that neither the basis 
set size nor the solvent drastically affect the spectral intensities. The 
computed spectral lines have been convoluted with Lorentzian 
functions of half-width 20 cm^^ See the main text for details. 


Gibbs free energy surface ^(Q). The calculation of if(Q) is 
nontrivial (cf., e. g., Refs. ZtII^ and citations therein). Nec¬ 
essary (but not sufficient) quantities that should be determined 
from electronic structure calculations are the adiabatic poten¬ 
tial energy surfaces En,a(Q) of the neutral and anion elec¬ 
tronic ground states at arbitrary geometries Q, as has been 
done in the present work. 


Applied to the case considered here, intramolecular reorga¬ 
nization manifests itself in the different locations Q ,4 ^ Q/v 


of the minima of the adiabatic energy surfaces Ea(Q) and 
Ea/(Q), respectively. The reorganization energies of the 
neutral and anionic species represent important quantities for 
electron transport both in the nonadiabatic (Mulliken-Hush- 
Marcus-type)^— and adiabatio^i^Likdl limits. They are de¬ 
fined as 


An = £'iv(QA)-£’ w(Qy), 
Aa = £a(Q;v)-£a(Qa). 


( 1 ) 


Our calculations confirm the intuitive expectation that sig¬ 
nificant contributions to the reorganization energy arise from 
the totally symmetric in-plane normal modes. The Raman ac¬ 
tive vibrational modes included in Table S3 belong to these 
modes. Out of them, the largest contribution to the reorganiza¬ 
tion energy comes from the mode 8a. Adiabatic energy curves 
En. A(G8fl) along the normal coordinate Q^a (all the other Q’s 
being set to zero) are presented in Fig. 0] and indicate a virtu¬ 
ally perfect harmonic behavior. 

Considered alone, each of the Raman active modes other 
than 8a listed in Table S3 has a small contribution to re¬ 
organization and also exhibits a harmonic behavior. Their 
contribution can be accounted for by a renormalization fac¬ 
tor of mode 8a //, (coga —>■ (Of, = fi, (Oga), as discussed in the 
Electronic Supplementary Information.— In accord with a 
general property of the harmonic modes,— the correspond¬ 
ing partial reorganization energies for the neutral and anionic 
species should be equal. For the mode 8a alone, we found 
= =0.0845 eV. 

However, the values of the total reorganization energies 
Ana computed from Eq.[T] were found different. An — 0.23 eV 
and Aa — 0.35 eV, and this inequality traces back to the low¬ 
est frequency mode, which is strongly anharmonic.— This 
mode is directly related to the inter-ring torsional motion 
and represents the floppy (label /) degree of freedom of 
the 44BPY molecule.— Its frequency cOf ~ 62 cm^*— lies 
well below the spectral vibrational range investigated ex¬ 
perimentally (400cm^^^a)^1800cm^') for 44BPY’’— and 
300 cm^^^a)~2000cm^*) for44BPY*^.~^ However, in spite 
of this very low frequency, its contribution to the reorganiza¬ 
tion energy is important because of the large amplitude vibra¬ 
tions associated with it. The anionic and neutral adiabatic en¬ 
ergy curves along the coordinate Qf of the neutral molecule, 
which are depicted in Fig. |5l exhibit strong anharmonicities. 
The fact that strong anharmonic effects characterize molecules 
possessing floppy degrees of freedom is well known (see, e. g.. 
Ref. l45h . However, except for the preliminary work of Ref. 
113 we are not aware of (transport-related) studies drawing at¬ 
tention to the fact that such modes can be strongly coupled 
to molecular orbitals with essential contribution to the charge 
transfer. 

In Fig.|6] we present adiabatic curves ENAiQf) for the neu¬ 
tral and anionic species calculated with various basis sets. The 
















Fig. 4 Adiabatic energy curves along the coordinate Qga of mode 8a 
of the neutral molecule (ro^ = 1640 cm^*) computed with the basis 
set aug-cc-pVDZ for neutral 44BPY and radical anion bipyridyl 
44BPY*^. Notice the almost perfect harmonic behavior and the fact 
that the neutral and anion minima are located on the same side with 
respect to the intersection point of the curves, situated at left beyond 
the figure’s frame (inverted regime). 



quencies presented in the Electronic Supplementary Informa¬ 
tion. On the contrary, the adiabatic curves of the radical anion 
exhibit a pronounced dependence on the basis set. A closer 
look at the -curves reveals that it is not the size of the 

basis set that matters; of primary importance is to include dif¬ 
fuse basis functions. For example, the results obtained by us¬ 
ing the cc-pVTZ basis set are much poorer than by using the 
rather modest 6-3lH-g(d) basis set. Noteworthy is also the fact 
that an anion unstable against electron autodetachment (i. e., 
Ea > Eff) can be the (incorrect) result of using an inadequate 
basis set; this aspect has already been pointed out.iii^ 



Fig. 6 The adiabatic energy curves £'a'a( 2/) the neutral 44BPYh 
and anion species along the coordinate 2/ of the floppy degree of 
freedom computed for several basis sets. For the neutral species, the 
curves computed with the basis sets 6-31-l-g(d) and aug-cc-pVTZ are 
practically indistinguishable. Notice the fact that, unlike £jv(2y), 
rather than the basis set size, the inclusion of diffuse basis functions 
is essential to correctly compute the adiabatic energy E^iQ/) of 
44BPY*^. Results obtained with the modest 6-31-l-g(d) basis set are 
reasonably good, while those with the basis set cc-pVTZ are not. 

See the main text for details. 


Fig. 5 Adiabatic energy curves along the coordinate Qf of the 
neutral molecule of the lowest frequency mode {(Of = 62 cm^*) 
computed with the basis set aug-cc-pVDZ for neutral 44BPY and 
radical anion bipyridyl 44BPY*^. Notice the strong departures from 
curves corresponding to the harmonic approximation. See the main 
text for details. 


inspection of these curves reveals an important difference be¬ 
tween the neutral and anionic species. The adiabatic curves 
of the neutral molecule are weakly affected by the basis set 
employed in the calculations; essentially, this dependence is 
similar to the geometrical parameters and the vibrational fre- 


5 44BPY in solvents 

We have also extensively investigated the adiabatic energy 
surfaces of the neutral (TTBEY**) and charged (radical an¬ 
ion 44BPY*^ and radical cation 44BPY*+) species in the 
solvents utilized in experiments, e. g., aqueous solution, 
toluene,— acetone,— 1,2,4-trichlorobenzene,— and acetoni¬ 
trile.— In line with the preliminary investigation of Ref. 
our calculations confirm the fact that, essentially, the solvent 
main effect is a nearly constant (i. e., Q-independent) shift 
with respect to the adiabatic surfaces of the corresponding 
species in vacuo. Therefore, we do not present adiabatic sur- 























faces in various solvents. Rather, we restrict ourselves to re¬ 
port results for the aforementioned energy shift, which is im¬ 
portant for molecular transport; the relative molecular orbital 
alignment with respect to the electrode Fermi levels is a key 
parameter that affects the current.— 

In Table [T] we present the computed influence of the sol¬ 
vents listed above on the vertical ionization potentials (IP= 
Ec — En = —Ehomo^ where the subscript C stands for radical 
cation) and the vertical electron affinities (EA= En — Ea = 
—Elumo) at the neutral equilibrium geometry. Although the 
44BPY-based molecular junctions fabricated so far do not 
seem to exhibit hole-type (HOMO-mediated) conduction, we 
have also included results for the radical cation 44BPY*+ and 
for the HOMO in order to point out an aspect, which is quite 
significant for the molecular transport in solvents. 

As visible in Table [T] the IP and EA-shifts are substantial 
and significantly dependent on the solvent. As a consequence, 
the HOMO-LUMO gap W/,; is significantly affected by the sol¬ 
vents. Instead of estimating this gap from the HOMO and 
LUMO Kohn-Sham (KS) “orbitals” (whose physical mean¬ 
ing is questionable^), in Table [T] we have determined it from 
A-DET calculations as W/,/ =IP—EA= Ec + Ea — 2En, an ex¬ 
pression which shows that it is the counterpart of what is called 
the c harge g ap in solid state and mesoscopic physics (see, e. g.. 
Refs. I 27 H 5 ll and citations therein). 

The point to which we want to draw attention here is that, 
although the solvent substantially affects the HOMO-LUMO 
gap, it causes HOMO and LUMO shifts in opposite directions, 
and their magnitudes are practically equal. The fact demon¬ 
strated here, that, for a given molecule (44BPY), a variety of 
solvents cause opposite HOMO- and LUMO-displacements 
of almost same magnitudes is complementary to another fact 
known earlier; the HOMO and LUMO displacements were 
found to be similar for a variety of (smaller) molecules in a 
given (namely, aqueous) solution.— 


6 Impact of the twisted molecular conforma¬ 
tion on the conductance 

As already noted in Ref. and reiterated above, unlike in 
the case of (practically) rigid molecules, the calculation of the 
electric current through molecular junctions based on floppy 
molecules requires a nontrivial step, namely an ensemble av¬ 
eraging in a nonequilibrium state.— This ensemble averaging 
is a problem on its own that will not be addressed in this paper. 

Instead, we will briefly consider a simpler problem related 
to it; the dependence of the ohmic conductance G on the 
twisting angle 6. Rather than examining changes in 6 due 
to molecular vibrations, in this section 9 will be considered as 
an independent geometrical parameter. This problem is also of 
interest; like in the case of the related biphenyl molecule, 


one can synthesize and investigate, e. g., derivatives of a 
molecular 44BPY series wherein the torsion angle 9 can be 
tuned by bridging the two pyridyl rings with alkyl chains of 
variable lengths. 

The dependence G = G{9), which has been considered 
theoretically for molecules consisting of two rings that can 
rotate relative to each other, interesting from the 

perspective of molecular conductance switching. Theoreti¬ 
cally, an approximate dependence G cos^ 9 has been found, 
which roughly agrees with experimental data in biphenyl- 
derivative This behavior characterizes a molecular 

switch that is “on” and “off” in parallel (0=0) and perpen¬ 
dicular (0 = 90°) conformations, respectively. 

In spite of significant differences between existing ap¬ 
proaches of the (ohmic) transport in 44BPY, which are 

based on nonequilibrium Green’s function (NEGE) methods 
and the Landauer trace formula (see, e. g.. Refs. 6ll-^), there 
is an important common finding in those studies. Namely, the 
fact that the coherent electron transport is determined by the 
LUMO: close to the metallic Eermi energy Ep, the calculated 
transmission exhibits a pronounced Lorentzian peak, which is 
well separated energetically from other peaks situated away 
from Ep. This single-level (“only LUMO”) model is simple 
but turned out to be realistic for 44BPY. By interpreting ex¬ 
perimental thermopower data within this model, the LUMO 
was found to lie by ^ 1.53 eV above the metallic Eermi en¬ 
ergy, in good agreement with the theoretically estimated value 
of ~ 1.47 eV.— So, we have a strong reason to accept this 
model, and then the low-bias conductance G can be expressed 
as (see, e. g.. Refs. 


.(e). 


Go £Ej/Mo(e)+r2 




( 2 ) 


Here E is a finite width due to the LUMO-electrode cou¬ 
plings,— which is usually much smaller than the LUMO- 
energy offset EpuMO relative to the Eermi energy Ep of the 
electrodes to which the molecule is linked. The quantity 
Elumo differs from the difference between the value Elumo of 
the isolated molecule and Ep, because of the energy shift 5E 
caused by molecule-electrode interactions (e. g., image charge 
effects^). We do not attempt to microscopically estimate 
here the difference between Elumo Elumo {Elumo = 
Elumo — Ep — 8E). We determine this shift by requiring that 
the value of EpuMoi^eq) at the equilibrium value 0 = 9eq 
(cf. Table SI) coincides with that deduced experimentally 
EpuMoi^eq) — 1-5eV for 44BPY junctions.— Based on the 
previous calculations,— which indicated that the metal (Au) 
atoms contacted to 44BPY essentially causes a constant {Qf- 
independent) shift of the LUMO energy (compare Eigs. 3 a 
and 3c of Ref. Wh . we will ignore the 0-dependence of 
Elumo {d) — EpuMoiS)- For the basis set aug-ccpVTZ, we 
thus deduced.Ei(7Mo(0) —EpuMoi^) — l-92eV. 

















Medium 

Whi 

IP= —Ehomo 

EA= —Eiumo 

SIP 

SEA 

SIP+S EA 

vacuum 

water 

toluene 

acetone 

acetonitrile 

TCB 

8.76; 8.69; 8.71 
5.26; 5.21; 5.24 
6.73; 6.67; 6.68 
5.40; 5.34; 5.37 
5.32; 5.26; 5.29 
6.82; 6.75: 6.77 

9.15; 9.13; 9.13 
7.42; 7.41; 7.42 
8.13; 8.12; 8.12 
7.48; 7.48; 7.49 
7.44; 7.44; 7.45 
8.18; 8.17; 8.16 

0.38; 0.44; 0.42 
2.15; 2.21; 2.19 
1.40; 1.46; 1.44 
2.08; 2.14; 2.12 
2.12; 2.18; 2.16 
1.35; 1.41; 1.39 

-1.73; -1.72; -1.71 
-1.01; -1.01; -1.01 
-1.66; -1.66; -1.64 
-1.70;-1.69;-1.68 
-0.97; -0.97; -0.97 

1.77; 1.76; 1.77 
1.02; 1.01; 1.01 
1.70; 1.70; 1.70 
1.74; 1.74; 1.74 
0.97; 0.97; 0.97 

0.04; 0.04; 0.06 

0.001; 0.0005; 0.003 
0.04; 0.04; 0.05 

0.04; 0.04; 0.06 

0.0001; -0.0006; 0.002 


Table 1 Values for the HOMO-LUMO gaps W/,; = Elumo^Ehomo^ ionization potentials IP, and electron affinities of the 44BPY molecule in 
vacuo and several solvents (TCB= 1,2,4-trichlorobenzene) computed with the basis sets 6-31-l-g(d); aug-cc-pVDZ; aug-cc-pVTZ, respectively. 
Differences between the values of a certain quantity in solvent and in vacuo are denoted by 8. Notice the opposite signs and the almost equal 
magnitudes 5 IP ~ —SEA of the changes in the ionization potentials and electron affinities caused by the solvents with respect to the values in 
vacuo. See the main text for details. 


Further, in view of the fact that changes of the molecule- 
electrode contacts turn out to only weakly affect the torsional 
angle,— we will also neglect a possible 0-dependence of F; 
we do not expect a substantial change in the short-range {con¬ 
tact) interaction between the LUMO and electrodes by vary¬ 
ing the inter-nt\g angle 0. By employing the experimental val¬ 
ues G(0eg)/Go = 0.00068 a.\\AEiijMo{Qeq) — 1.5eV,— Eq. (|2]) 
yields the value F ~ 0.039 eV. 

To obtain the LUMO energy Elumo{(^) of the isolated 
molecule, we performed A-DFT-calculations— and com- 
putedH 

ELc/Mo(0) = EA(0;R')-£iv(0;R')- (3) 


For simplicity, we present results of the calculations carried 
out for various conformations by varying the torsion angle (0) 
while freezing all the other geometrical parameters at the neu¬ 
tral equilibrium geometry (denoted by R').— Results illustrat¬ 
ing the 0-dependence of the LUMO energy Eujmo{Q) of the 
isolated molecule are shown in Lig. [T] 

Once Eujmo{Q) is known, the 0-dependence of the con¬ 
ductance can be deduced by means of Eq. (|2]i. The depen¬ 
dence of G on cos^ 0 is shown in Lig.|8] As visible there, our 
result does not give support to the proportionality G cos^ 0 
found experimentally for other molecule s and DLT-based 
theories approaches. 35,58,60 q^j- results agree very well with 
the following Ansatz, which straightforward generalizes the 
aforementioned relationship 

G(0)/Go=gx+(g||-gx)cos2 0. (4) 


It is worth emphasizing at this point the fact that in 44BPY 
the relationship G cos^ 0, which yields a perfect “off” 
switch G(0 = 90°) = 0, is contradicted by a detailed micro¬ 
scopical analysis.— According to Ref. the transmission 
through 44BPY remains nonvanishing even in a strictly per¬ 
pendicular conformation, due to the existence of crossed a — K 
electronic interaction.— According to the present estimation, 

^±-^||/ 3 - 





Fig. 7 A-DFT results for the LUMO energy of the 44BPY molecule 
computed with several basis sets (specified in the inset). Notice the 
importance of the diffuse functions. See the main text for details. 


We end this sections with two remarks. Lirst, we empha¬ 
size the difference between our present A-DLT-based calcu¬ 
lations and most DLT-approaches, wherein mathematical ob¬ 
jects, namely KS orbitals, are treated as if they were real 
physical orbitals. Second, we note the relationship of the 0- 
dependence discussed above and the very small energy bar¬ 
rier associated to the profile of the torsional potential, which 
is Eiq{9) in our notation. An important issue amply dis¬ 
cussed in the past (see, e. g.. Refs. 18lfl9 and citations therein) 
is the magnitude of the barrier heights between the twisted 
global minimum and the planar and perpendicular conform- 
ers, AEq = En{0 = 0°) -£A,(0e^) and AEgo = Eig{9 = 90°) - 
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Fig. 8 Dependence of the conductance G in units of conductance 
quantum Gq = 2e^/h = 77.5 on the torsional angle 6 obtained 
hy using the basis sets specified in the inset. For each basis set, the 
conductance has been been fixed at the experimental value 
G{eeg)/Go = 0.00068ii (Beg = 37.01 for aug-cc-pVDZ and 
Beg = 37.3 for aug-cc-pVTZ, cf. Table SI). The straight lines have 
been obtained by fitting the results of quantum chemical 
calculations using the Ansatz of Eq. 0. The fact that the symbols 
corresponding to the various basis sets employed are hardly 
distinguishable from each other demonstrates the realiability of the 
results. See the main text for details. 


EN{Seg), respectively. The most recent calculations published 
so far— succeded to produce values A£o —1.5 kcal/mol and 
AEgo ~ 2.2 kcal/mol, which agree with the ordering AEq < Ago 
deduced from measured NMR spectra of 44BPY dissolved in 
nematic liquid crystal.— However, there still remain signifi¬ 
cant differences between theory and experiment; proton NMR 
measurements estimated AEq ~ 4.0kcal/mol.^' Internal rota¬ 
tional barriers of heights comparable to 44BPY have also been 
found for other molecular systems. 


7 Discussion and Conclusions 

In this paper, we have reported results of extensive quantum 
chemical calculations for the 44BPY molecule in vacuo and 
in different solvents employed in relevant experiments. 

The small differences between the measured properties of 
the neutral 44BPY‘' and radical anion 44BPY*^ and those, 
nearly convergent with respect to the basis set size, calculated 
within the DFT/B3LYP presented in the Electronic Supple¬ 


mentary Information are representative for the present state- 
of-the-art of the theory: on one side, it demonstrate its rather 
high accuracy and, on the other side, that intrinsic limitations 
remain. 

The results on the adiabatic energy curves and on the sol¬ 
vent effects presented in this work can (and will) serve as input 
information for subsequent transport studies through 44BPY- 
based molecular junctions. 

Because, as expressed by the title, the remote goal of the 
present investigation is the molecular transport, three findings 
of this work are particularly worth to emphasize. 

(i) We are not aware of attempts to microscopically incor¬ 

porate solvent effects into transport approaches. Available 
program packages do not allow to perform microscopic trans¬ 
port calculations through molecular devices immersed in elec¬ 
trolytes. Our results indicate that the solvent acts as a selec¬ 
tive gate electrode: it causes energy shifts of particle-like and 
hole-like orbitals in opposite directions and nearly same mag¬ 
nitudes. On this basis, we propose to embody solvent effects 
on molecular transport within a procedure based on a scissor 
operator. Let us remind at this point that the scissor operator 
technique has been proposed— and applied—"— in semicon¬ 
ductors to empirically cut the band structure along the band 
gap and rigidly move the energies of the conduction and va¬ 
lence bands in opposite directions. There, the scissor operator 
has been motivated technically, as an (empirical) correction 
of the band gap, whose value is underestimated by the DFT. 
More recently, a scissor operator has been utilized in molecu¬ 
lar electronic s with both a technical and a physical mo¬ 

tivation: to correct the too narrow KS HOMO-LUMO gap and 
to account for the energy shift due to electrodes via image ef¬ 
fects,— respectively. For simplicity, in those approaches the 
HOMO and all occupied orbitals, and the LUMO and all un¬ 
occupied orbitals are shifted in opposite directions by the same 
quantity A/2 (half the correction to the HOMO-LUMO gap). 

Similar to those discussed above, we also propose to ac¬ 
count for solvent effects on the molecular transport by means 
of a scissor operator: the energies of the occupied and un¬ 
occupied orbitals should be moved by —5 IP and by —5 FA, 
respectively. Unlike in the aforementioned cases, wherein the 
assumption that the opposite shifts have the same magnitude 
(A/2) is merely done for convenience, our results give a mi¬ 
croscopic support to the fact that the shifts are opposite and 
have practically the same magnitude —5IP~ -t-5EA. Notice 
that (a) 5 IP and 5 FA are not empirical quantities but deduced 
from quantum chemical calculations, and (b) the presently 
proposed scissor operator procedure can be applied on top of 
any transport approach in vacuo, let it be a ‘DFTh-E” imple- 
mentatio n or else. 

(ii) The results presented here (Figs. |6] and |7]i clearly 
demonstrate that, in order to correctly describe the anion 
(thence the LUMO), it is necessary to employ basis sets that 








include diffuse basis functions. Even by using basis sets of 
double zeta quality (cc-pVDZ), the omission of the diffuse 
functions, as often the case in many transport calculations, 
may yield errors ^ 0.4 — 0.5 eV in the LUMO energy off¬ 
set. In view of the fact that typical energy offsets amount 
^ 0.6 —1.5 eV, such errors in determining the LUMO 

location are unacceptably large. 

(iii) Along with molecular junctions wherein a single 
44BPY molecule is directly linked to electrodes, 
molecular junctions wherein spacers (e. g., six alkyl units at 
either side-i^) are used to link the 44BPY molecule to elec¬ 
trodes are also of experimental interest. Because the largest 
basis sets employed in the present work may become a chal¬ 
lenge for such molecular sizes, it is important to note that, 
once diffuse basis functions are included, anionic species can 
still be reliably described with rather modest basis sets, as il¬ 
lustrated by the results obtained by employing the 6-3lH-g(d) 
basis set shown in Figs. |6] and |7] and Table[i| 

Figs.|6]and|7]present results demonstrating both the impor¬ 
tance of (a) the diffuse basis functions and of (b) the most 
salient structural feature of the 44BPY molecule, namely the 
twisting inter-ring angle. One should emphasize that these two 
aspects are distinct from each other. Concerning aspect (a), 
one should note that diffuse functions are needed for the cor¬ 
rect description of the anion. The errors introduced by using 
basis sets without diffuse functions, e. g., at the neutral equi¬ 
librium geometry = 0 in Fig.|6]or 0 = 9eq — 37° in Fig.|2l) 
are practically the same as those far away from it (cf. Figs. |6] 
and|7J. Concerning aspect (b), we remind what was already 
stated in Ref. (see especially the next-to-last section of that 
work): the floppy degree of freedom (Qf) as well as its strong 
anharmonicity and significant reorganization energy are prop¬ 
erties of the 44BPY molecule, which are affected neither by 
solvents nor by electrodes. Computations for 44BPY'^ and 
44BPY*^ immersed in solvents and by attaching gold atoms 
at the two ends show that the main effect is a practically con¬ 
stant shift of the anionic adiabatic energy surface with respect 
to the anion placed in vacuo. 

To account for the effect of the intramolecular reorganiza¬ 
tion related to the floppy degree of freedom on the electron 
transport through 44BPY-based junctions, one should perform 
an ensemble averaging on a molecular system at away from 
equilibrium (nonvanishing external bias), which is a nontrivial 
problem.— As a preliminary step in demonstrating the signif¬ 
icant impact of molecular conformation on the transport, we 
have presented results for the linear (ohmic) conductance as 
a function of the torsional angle G = G{9). Our results ob¬ 
tained by quantum chemical calculations are well described 
by a phenomenological Ansatz, Eq. ©. 

To end, let us finally mention that, at 0 = 90°, the 44BPY 
molecule (02^; symmetry) becomes unstable against a Jahn- 
Teller distortion. Therefore, significant deviations from 


the above Ansatz may occur close to the perpendicular con¬ 
formation. Obviously, this is an issue that escapes the present 
considerations. 
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